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Abstract 

The recent proposal connecting the AdS/CFT correspondence and entanglement renormalization 
tensor network states (MERA) is investigated by showing that the entanglement entropy and the 
two point functions in a type of hybrid tensor network state composed by a finite number of 
MERA layers and a matrix product state (MPS) acting layer, imitate the behaviour of 

the holographic entanglement entropy and the two point functions in the AdS soliton geometry. 
Within the context of AdS/CFT, AdS solitons represent theories with a mass gap, i.e gapped 
systems. From these observations, an explicit connection between the entanglement structure of 
the tensor network and those parameters which define the AdS soliton geometry is provided. 
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I. INTRODUCTION 

The best known concrete realization of the holographic principle is the AdS/CFT corre- 
spondence l|. This duality is generally formulated as the equivalence between the generating 
functionals Z of two very different theories: a classical theory of supergravity defined on 
a D + 1 asymptotically Anti de Sitter spacetime and a QFT living on a D-dimensional 
spacetime that is the boundary of the AdS, i.e, 

Z QG [AdS D+1 , J(x, z)\ = Z QFT [d(AdS) D , J(x)} , (1) 

and more explicitly, 

e -S QG [J(x,z)} ^ / -JJ(x)0(x)d D x\ / 2 N 
\ / J(x,0)=.J{x) ' 



2 



where x are the spacetime dimensions of the QFT, z is the extra-dimension of the gravity 
theory on AdS, O are operators of the quantum field theory on the boundary, and 3{x, z) 
are the classical fields of the gravitational theory whose boundary values J(x, z)\ z ^ act as 
the sources J(x) for the correlation functions of the QFT. 

Since the initial formulations of the correspondence, a huge amount of work has been 
carried out in order to generalize the AdS/CFT while providing numerous valuable models 
to study non-perturbative effects in quantum field theory such as confinement and quantum 
phase transitions. However, still there is not a first principles derivation of the correspon- 
dence that allow us to explicitly construct a bulk gravity theory from a boundary QFT . 
Nevertheless, it is widely accepted that the AdS/CFT is at heart a geometric formulation 
of the renormalization group (RG), in which the renormalization scale becomes the extra 
radial dimension z and the beta functions of the boundary field theory are the saddle point 
equations of motion of the bulk gravity theory [3J. 

In recent years, ideas coming from quantum information have been especially relevant 
in establishing new numerical real space quantum renormalization group methods such as 



density matrix renormalization group (DMRG) Q] 
rithms including matrix products states (MPS) 



and tensor network states (TNS) algo- 



m, projected entangled-pair states (PEPS) 
6], multi-scale renormalization ansatz (MERA) 9j, tensor renormalization group (TRG) 



7|, and tensor-entanglement-filtering renormalization (TEFR) [8j]. These algorithms pro- 
vide a set of variational ansatze useful to characterize the low-energy (long distance) physics 
of quantum many-body systems. The TNS assumes a parameterization of a many-body 
wave-function by means of a collection of tensors connected into a network. The number of 
parameters required to specify these tensors is much smaller than the exponentially large 
dimension of the system's Hilbert space, allowing for the efficient representation of very large 



(and even infinite) systems. It has been recently proposed 11( that tensor network states 
can be broadly classified into two categories according to the geometry of the underlying 
networks. 

In the first category, the network mimics the physical geometry of the system, as specified 
by the pattern of interactions in the Hamiltonian. The MPS ansatz for one dimensional 
systems, which consists of a collection of tensors connected into a chain and its generalization 
for higher dimensional systems (PEPS), lie on this physical geometry category. 

At variance, according to ll|, the collection of tensors in the second category of TNS 
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are connected so as to parametrize the diferent length scales (or, equivalently, energy scales) 
relevant to the description of the many-body wave-function. These tensor networks organize 
the quantum information contained in a state in terms of different scales through a charac- 
teristic tensor connectivity which spans an additional dimension related with the RG scale, 
as for instance in MERA. This has been used to define a generalized notion of holography 
inspired by the AdS/CFT duality [10]. There, the author has made the groundbreaking ob- 
servation that the tensor networks in MERA, happens to be a realization of the AdS/CFT 
correspondence. In |loj|, the dual higher dimensional holographic geometry of the original 
system emerges when one realizes that the tensors in MERA corresponding to a one di- 
mensional quantum critical point, are connected so as to reproduce a discrete version of the 



anti-de Sitter spacetime (AdS). After the initial proposa 



a substantial amount of work has 



I1M22]. 



appeared supporting and extending the original idea 

Within the context of the AdS/CFT correspondence, the holographic dual of a gapped 
system is the capped holographic geometry called the AdS soliton This capping of the AdS 
geometry implies an infrared (IR) fixed point at a finite energy scale for the original CFT 
lying at the boundary of AdS, thus yielding the dual for a massive deformation of the theory. 
As the AdS soliton is dual to a gapped system, it is expected that the IR state will be a 
non entangled state if there is no topological order. However, it also might be the case of 
nontrivially entan gled IR states in extremal solitons, which could be due to some underlying 
topological order [16] . In [l6j, authors speculated, in terms of the AdS/MERA duality, 
about the structure of a tensor network that might suitably describe this possibility. 

The aim of this paper is to offer new insights on the connection between the structure of 
MERA states and their potential holographic descriptions by exploring some hints posed in 
1g| . Here we explicitly show that the entanglement entropy and the two point functions in 
a type of hybrid tensor network state (composed by a finite number of MERA layers and 
an MPS acting layer) representing a one dimensional quantum many body system, 

imitate the behaviour of the holographic entanglement entropy and the two point functions 
in the capped holographic geometries described by AdS solitons. From these observations, 
we provide a connection between some parameters of the tensor network representation of 
the system and those parameters defining the holographic geometry of an AdS soliton. 

The paper is organized as follows: In section two, MPS and MERA tensor network 
states are briefly reviewed and the hybrid tensor network composed of a finite number of 



4 



MERA layers and a capping MPS on the top layer is presented. Section three reviews, 
the AdS/MERA conjecture and its relation with the formula for the computation of holo- 
graphic entanglement entropy [23 ] . In section four, the entanglement entropy and the two 
point functions in AdS soliton geometries are contrasted with those results obtained when 
computing these quantities in the hybrid tensor network presented in section two. Finally, 
in section five, results are summarized and future research problems are exposed. 



II. TENSOR NETWORK STATES FOR 1-D SYSTEMS 

Tensor network states constitute a new class of numerical methods that aim to efficiently 
describe ground states of strongly correlated quantum systems. TNS implement techniques 
drawed from renormalization group methods and the knowledge about the structure of 
entanglement in the ground states of quantum many body systems. In this section we 
briefly review some issues concerning the two classes of tensor network states (MPS and 
MERA) used in this paper in order to offer further evidence of the connection between 
AdS/CFT and MERA. 



A. Matrix Product States 



Let us first to consider a 1-D translationally invariant quantum many body system com- 
posed by N sites of dimension d. Its ground state can be expressed as 

d 

1^) = ^ %i,s 2 ,~ ,*n\ S 1> S 2, • ■■ ,S N ) , (3) 

si,S2,— ,sjv =1 

where for j G {1, - • • , N} and Sj G {1, • • • , d}, the vectors \sj) G Hd form the com- 
putational basis of the j-th system site and where the type - (^j tensor % lt s 2 ,-,s N = 
(si,S2, ■■■ jSjvI^) represents the associated probability amplitudes. The MPS representa- 
tion of this state is given by the ansatz {5], 

T Sl , S2 ,-,s N — Tr (A Sl A Ss ■ ■ -A SN ) , (4) 

where A Sj are matrices of dimension M. with elements A S J ^ whose values are determined 
by means of a variational search procedure. Away from the critical point, a finite A4 suffices 
to exactly describe the ground state. Essentially, the so called bond dimension A4, controls 
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FIG. 1: Pictorial representation of the MPS ansatz. The matrices A Sk are three indices tensors 
(excepting the first A Sl and last A SN tensors, which only have two indices due to open boundary 
conditions), with Sk related with the physical site k and a and (3 as the indices of the correlation 
space of dimension M linking to the two neighboring matrices A Sk ~ 1 and A Sk+1 . In the figure, 
connecting two adjacent tensors through a given link is equivalent to contracting the product of 
these two tensors through the corresponding index. Applying this rule to all the matrices of the 
picture results in equation (TIJ). 

how much entanglement between adjacent sites is retained during the variational procedure. 
Treating the matrices A S3 as f * J tensors (one extra index is due to the physical site), the 
MPS ansatz can be pictorially described as in figure 1. 

Regarding the computation of observables with the MPS ansatz, given the state ([3]) with 
coefficients (J4]), the expectation value of an operator = 0\ ® 2 <8> • • • ® On which is the 
tensor product of local operators Oj for each site j, can be expressed as (Fig. 2), 

(M>|0|M>) = Tr(E Ol E O2 ■••E ON ) . (5) 

The transfer matrices K . are denned by, 

d 

E o 3 = £ (s'jio^) (,k< ,r') , (6) 

and A s j is the hermitian conjugate of A Sj . As a result, as posed in equation (jSJ), the 
computation of expectation values with MPS, may be viewed as the total contraction of a 
1-dimensional array of tensors (Ej). 

It is well known that an MPS with a finite bond dimension Ai, faithfully represents the 
ground state of a ID gapped system and supports the exponential decay of correlations 
expected for these models. Considering a translational invariant matrix product state |\&), 
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FIG. 2: Pictorial representation of an expectation value computation through the MPS ansatz 
given by equations © and ([6]). 

homogeneously denned by the matrices A s , the two point correlation function of two local 
observables Q and 0/3, with support on sites Sj and Sj separated an arbitrary distance £ + 1 
is defined by, 



<£g 1] = (*|6 a ( S< ) 0^)1*) = Tr (Ef l] • • -Eg • • -Eg • • -E? 



Following 24j . the correlator in the limit when iV — >■ oo can be written as: 



r [m] 

'a/3 



E 

i/=2 



(7) 



(8) 



where A„>2 are the eigenvalues of E x for which it holds that | A^>2 1 < 1, \v) and (y\ are the 
right and left eigenvectors of E x for the eigenvalue A„, V^ v \t) = (1|Eq^|z/) (z/|Eq^|1) are 
polynomial functions of degree H u with coefficients depending on X u , and B v is the size of 
the largest Jordan block belonging to the generalized eigenspace of A„. 

Since every |A^> 2 | < 1, then the correlator V^t decay exponentially because it is possible 
to write its leading behaviour as a superposition of exponentials with decay lengths defined 
by £„ = — l/log|A w |. This implies that the MPS gapped nature dominates when probing 
long range correlations ruled by the eigenvalues \ u of the identity transfer matrix Ej. 

Furthermore, if one computes the entanglement entropy of a subsystem A within an 1D- 
MPS state, bipartitioning the system into two subsystems A and B, and then taking the 
partial trace over the sites on B, it is easy to convince oneself that, 



S A < 2 losM. 



(9) 



Thus, Sa does not scale with the size of the subsystem A despite it can be made arbi- 
trarily large by increasing Ai. As the entanglement entropy in ID critical systems scales 
proportional to log£, where £ is the size of the subsystem A [39, 40], then, from the bound 
one clearly sees that an MPS does not support the large amount of entanglement required 
in quantum critical systems. 

B. Entanglement renormalization tensor network states. MERA 

The MERA representation of the state l-D state §3§ assumes a decomposition of 
T ai ,82,-,3N * n terms of a collection of smaller, finite size tensors, which differently from the 
linear MPS structure, are organized in a two-dimensional layered graph. The sites of the 
graph represent tensors. They are divided in two groups: the type \ ~ ( \ J tensors with 
elements Xsi,'" 2 2 caue d disentanglers and the type A — ( * J tensors with elements A" i 1 S2 called 
isometries (Fig. 3). For quantum critical systems, a MERA tensor network posses a charac- 
teristic scale invariant structure, i.e, a unique \ and a unique A precisely define the MERA 
graph. 

The MERA representation of ([3]) implements an efficient real space renormalization group 
procedure through a tensor network organized in different layers labelled by r. Each layer 
of MERA defines a RG transformation: prior to the renormalization of a block of typically 
two sites located at layer r into a single site by means of a A T -type tensor, short range 
entanglement between the sites is removed by means of the disentangler \ T ■ In other words, 
the MERA tensor networks, as a consequence of the locality of physical interactions in the 
ground state of many-body systems, at each layer, decouple the relevant low energy degrees 
of freedom from the high energy ones, which are then safely removed, by unitarily transform- 
ing with disentanglers small regions of space. The MERA coarse-graining transformation 
induces an RG map that can be applied arbitrarily many times keeping constant the com- 
putational cost for obtaining consecutive effective theories for the original system which are 
also labelled by r. 

In a ID MERA the computation of an observable G requires the contraction of a 2D 
tensor network, and this may be efficiently done by the isometricity requirements for \ and 
A tensors (see figure 4). Indeed, in an scale invariant MERA, contracting an entire layer of 
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FIG. 3: Pictorial representation of two layers (r and r + 1) of an scale invariant MERA tensor net- 
work. Three legged triangles represent isometries A and four-legged circles represent disentanglers 
X- At each layer r, the tensors are chosen so as to fulfill x X = I an d A' A = I. 



MERA and its conjugate, maps the original expectation value problem into a new one i.e, 

(*|e|*) = |tf'> = (^'|e [1] |^') , (io) 

where the quantum state |\&') is the original MERA state \^f) after the full contraction of the 
tensors at the bottom layer with their conjugates, the observable G)W = S (6) is the effective 



tensor shown inside the rectangle in figure 4 (right), and S is the scaling superoperator [41] 
of the observable B. The total contraction of the tensor network proceeds by mapping the 
effective observable in a sequential way as, 

qM e [r+i] = 5(e M) ? 

until the top of the tensor network \C) is reached, for which it holds that (\l/|0|\l/) = 
(C|0^|C), with h = log 2 N and being the total size of the system. 

One of the most salient features of an scale invariant MERA tensor network is that it 
naturally supports the well known power-law decaying of correlation functions in systems at 
the quantum critical point [26]. To show this, it is convenient to define the fixed point one 
site scaling superoperator S = as, 

5(i ) (o) = ^ /UQ $ a Tr($ Q o) , (12) 

a 

where the scaling dimensions of the scaling operators $ a of the theory are A a = — log 2 [i a 
and Tr($ Q $#) = 5 a p. Let us consider the correlator of two scaling operators <& a (sf) and 
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0-1 

FIG. 4: To compute the expectation value of an observable (blue squares) supportted on an 
interval of sites, {si-.s*-} with a ID MERA state, the first layer of disentanglers f° r which, 
by definition, holds x X = I is initially considered. Those Xn'r™ 2 whose (both) physical indices do 
not belong to the support interval {si. are trivially contracted with their adjoints. However, 
this contraction is not trivial for those disentanglers directly connected to the support. The same 
argument holds for the isometries A belonging to the first layer of the MERA state. 

®p(sj) initially located at sites Si and Sj of the original lattice, 

e^H*l*«(*)**fo)l*> . (is) 

After r* = log 2 \s{ — Sj\ contractions of the type described above, the sites supporting the 
original $ a and become first neighbors [26| and each iteration contributes a factor ji a fip 
giving, 

(H?\* a {8 i )9p{ 8j )\9)= (14) 

= (c*\(s {l) ^ a ( Si )Y {s w (Msi)Y in = 

= (/i Q ^) T (C*|$ Q (0)^(l)|C*) , 

where the state \C*) is the result of r* contractions of the original MERA state with 
its conjugate. Defining £^ = (C*|$ a (0) <3>g(l)|C*}, noticing that A aj/3 = — \og 2 n a ,p an d 
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recalling the identity a g = b ga , the two point function (TT3I) can be written as, 

^ a ( Si )^( Sj )) = -^-, (15) 

I I 

with 77 = (A Q + Ajg)^ = 2A Q . 

C. Hybrid MERA-MPS networks 

When considering the ground state |\?) of gapped lD-Hamiltonian "H, the correlations 
decay exponentially at large distances due to the characteristic correlation length £, while 
typically keep power-law decaying for distances smaller than £. To emulate this behavior 
with a tensor network state, it is necessary to have an ansatz which suitably reproduces these 
features. A potential candidate for this, is a tensor network state with an MPS at the top 



of a finite number (tq ~ log£) of MERA layers U|. Here it is assumed that these To layers 
correspond to those in the scale invariant MERA describing a neighbouring critical point 
while the capping MPS accurately describes the behaviour of the long distance correlations 
in the gapped phase. 

In the MERA representation of the ground state of a gapped Hamiltonian, after r ~ log £ 
renormalization steps, the original ground state |\I/o) of the system has flowed into a state 
\Q) that can be well approximated by a product state with no entanglement between the 
different coarse grained lattice sites. This state describes the ground state at a fixed-point 
of the RG flow corresponding to a gapped phase without topological order. Nevertheless, 
the finite layered MERA may also be combined with another tensor network (an MPS in 
figure 5) in order to represent the ground states of gapped systems that flow towards non 
trivially entangled IR fixed point ground states |$). 

The MPS tensor at the top layer of this hybrid MPS-MERA state, accounts for the 
gapped character and non trivial long range entanglement of the IR fixed point in a natural 
way, i.e, it naturally arranges the exponential decay of the two point correlations at distances 
larger than £, while the MERA 'curtain'' of To layers implements the power-law decaying 
behaviour of correlators when the distance between the physical sites is smaller than £. The 
case of the trivial non-entangled product IR state \ is represented by an MPS with bond 
dimension M. = 1. 
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FIG. 5: Pictorial description of an Hybrid MPS-MERA tensor network state with a fixed number 
To of MERA layers and a top tensor given by an MPS. The original ground state of the system 
|^o) nes & t the first layer of the MERA curtain while the ground state at the IR fixed point is 
given by the MPS tensor capping the MERA geometry. 

In the MPS-MERA tensor network, the total system size N is given by, 



where m, is the number of effective MPS sites located at the top of the To-th layer of the 
binary MERA ' curtain 1 hanging from the MPS tensor. Regarding the entanglement that 
may be supported by these type of tensor states, as each "site" in the top MPS represents 
a cluster of ~ 2 T ° coarse grained sites of the original lattice, then the entanglement between 
different sites of this tensor represents the long range entanglement between distant clusters 
of sites in the original system. As a result, if the MPS cap tensor is not the trivial product 
state, its long-range entanglement, which relates with some type of topological order, will 
have a non zero contribution to the total entanglement entropy of a region A in the original 
lattice. As an example, the ground states of ID systems in the Haldane phase, show a kind 
of symmetry protected topological order [271 ] ■ In the tensor network representation of a 
system with this type of topological order, the local unitary transformations carried out by 
the disentanglers of the MERA "curtain" only remove the short range entanglement, i.e., 
the entanglement between the neighboring sites. At variance, the long range entanglement 
(i.e, the entanglement between distant sites), is not removed by disentanglers. It is precisely 
this robustness of the long range entanglement under local operations which entails the non 



N 



rn 



(16) 
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trivial topological order of the IR fixed point state represented by the MPS at the top of 
the MERA curtain. 



III. THE AdS/MERA DUALITY 



19 



In [10||, it was firstly observed that MERA [9[ happens to be a realization of the AdS/CFT 
correspondence 1]. As pointed out in [lo|, from the entanglement structure of a quantum 
critical many body system, is possible to define a higher dimensional geometry in which, 
apart from the coordinates labelling the position and the time t, one may add a "radial" co- 
ordinate z labelling the hierarchy of scales. Then, the higher dimensional geometry emerging 
from MERA may be visualized by locating cells around all the sites of the tensor network 
representing the quantum state. The size of each cell is defined to be proportional to the 
entanglement entropy of the site in the cell. Through this procedure, a geometric dual pic- 
ture of the tensor network arises quite naturally from the entanglement of the degrees of 



freedom of the critical system lying on the boundary [12|, |20 |. 

The discrete geometry emerging at the critical point is a discrete version of AdS. For a one- 
dimensional quantum critical system with a space coordinate labelled by x, the continuous 
isometry w — > w + 9, x — > e e x of the metric 

ds 2 AdS ~ dw 2 + e~ 2w (dx 2 - dt 2 ) , (17) 

is replaced by the MERA's discretized version, w — > w+k, x — > 2 k x or x — > 3 k x depending 



on the binary or ternary implementation of the renormalization algorithm [9]. Here, the 
extra direction w in the AdS is identified with r (the variable labelling the number of 
renormalization steps) in MERA. To be more concrete, the (D + l)-dimensional AdS metric 
can be written in the form, 

ds l dS = Sf. ( d z 2 - dt 2 + dx 2 ) , (18) 

where, LAds is a constant called the AdS radius; it has the dimension of a length and it 
is related with the curvature of the AdS space. We have also made the change of variable 
t = logz, where z is the standard radial coordinate of the Poincare AdS. With this choice of 
the spacetime coordinates, the one dimensional quantum critical system lies at the boundary 
{z — 0) of the bulk geometry. 
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The most convenient way to establish the AdS/MERA connection is to compare the 
procedures to compute the entanglement entropy in both cases. In the classical gravity 
limit of AdS/CFT, Ryu and Takayanagi (RT) derived a celebrated formula yielding the 
entanglement entropy of a region A provided that the (boundary) conformal field theory 



describing the critical system admits an holographic gravity dual |23J]. In the RT approach, 
the entanglement entropy is obtained from the computation of a minimal surface in the 
dual higher dimensional gravitational geometry (bulk theory); as a result, the entanglement 
entropy Sa in a CFT/? is given by the area law relation, 

_ Areafa) 

SA ~ 13™ ' ( 5 

where D is the number of spacetime dimensions of the boundary CFT, 7^4 is the D- 
dimensional static minimal surface in AdS^+i whose area is given by Area(7^) and G^ +1 ^ 
is the D + 1 dimensional Newton constant. In the RT proposal, looking for the minimal 
surface 7a separating the degrees of freedom contained in region A from those contained 
in the complementary region B, amounts to search for the severest entropy bound on the 



information hidden in the AdS^j+i region related with B. Despite the RT formula 



ras not 
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28 



29]. 



been rigorously proven its validity is supported by very reassuring evidence 

In this paper, the interest will be mainly focused in the D — 2 = (1 + 1) case, for which 
the equation (Tl9l reduces to, 

= Length^) 
4G< V 3 > 

Namely, there is a striking similarity between the formula ( )20|) and the computation of the 
entanglement entropy of a region A with MERA. Using MERA, the entanglement entropy 



is fixed by the structure of the causal cone 



25]. The causal cone CC(A) of a region A 



of I sites, is determined by all the disentanglers and isometries along all the levels of the 
tensor network, which are directly connected with the sites in the original region A, (see 
figure 6). As a result, to compute the entropy Sa it is necessary to trace out any site within 
the tensor network which does not lie in CC(A). The boundary of the CC(A) is a curve 74. 
The length of 74 is, when multiplied by the maximal contribution to entanglement entropy 
given by each traced out site, the sum of the entropies of all the traced out sites 0, [hJ, and 
thus provides an upper bound for the entropy Sa [11]. 

S A < Length^) log W , (21) 
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FIG. 6: The causal cone of a block of sites in MERA (in red) is given by the minimal number 
of tensors that may affect to the expectation value of any observable with support on those sites. 
The boundary of the causal cone of an interval A defines a curve 74 which length, in an optimized 
MERA, acts as the bound for the entanglement entropy of the region A given in equation (|21|) . 



where W is the dimension of the auxiliary space of disentanglers and isometries. For 
instance, for disentanglers Xvlvh M i,2 £ {1| • " i ^} an d ^1,2 £ {1, • • • , W}. 

The close connection with the RT formula emerges as one realizes that 7^4 can be regarded 
as the minimal curve 7^ in the RT proposal, since it counts the minimal number of sites 
which must be traced out in the MERA tensor network. Indeed, the minimal curve m 



an optimized scale invariant MERA network llj saturates the bound given in (12 lj) : this has 



been confirmed by explicit computation in one dimensional critical systems, where it has 



been shown that 



9], 



S A ~k-\og£. (22) 



where A; is a constant of order one. Since 7,4 is defined as the boundary of the CC(A), it can 
be interpreted as an holographic screen which optimally separates the region in the MERA 
tensor network related with the degrees of freedom of A, from those in its complementary 
region B. 

Despite a complete characterization of the AdS/MERA duality is still lacking, very re- 
cent advances point to fulfill this objective. In j^j, authors have used a continuos version 



of MERA (cMERA, 2lj), the holographic formula for the entanglement entropy and the 
concept of distance between quantum states to address one of the major issues concerning 
the AdS/MERA conjecture, i.e, how gravity duals arise from the structure of the quantum 
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entanglement of the quantum critical systems under consideration. The question about the 
role of the large N limit in the AdS/MERA duality is also considered in (2^] as well as in 



221 ]. where it is shown how a number of features of the holographic duality in the large N 



limit emerge naturally from entanglement renormalization tensor networks. 



IV. HOLOGRAPHIC ENTANGLEMENT ENTROPY IN AdS SOLITONS AND 
HYBRID MERA NETWORKS 

It is widely accepted that, non-local observables such as two-point functions, entangle- 
ment entropy and Wilson loops turn out to be crucial for elucidating the subtleties of the 
AdS/CFT correspondence jlj. In this section, the holographic entanglement entropy and 
general features of two point correlators in AdS soliton geometries are calculated while inter- 
preting the results in the light of those obtained when computing these quantities by means 
of the hybrid MERA tensor network described in section III CI 

Within the context of the AdS/CFT correspondence, the holographic dual of a gapped 
system is the AdS soliton, obtained by double Wick rotating the AdS black hole and then 
compactifying one of the spatial dimensions transverse to the AdS radial dimension. The 
compact dimension shrinks to zero at some finite value zq of the AdS radial coordinate, 
capping the AdS geometry. The radial coordinate together with the shrinking of the compact 
dimension make a type of cigar geometry, with the IR end of space corresponding to the tip 
of the cigar. Furthermore, this capping of the AdS geometry implies an IR fixed point at 
a finite energy scale (~ l/z ) for the CFT lying at the boundary of AdS, thus yielding the 
dual for a gapped massive system. 

Here it will be considered the AdS3 soliton metric, which in Poincare coordinates is given 

by, 



ds 



T 2 

Sol" 2 



dz 



2 



+ f(z) dx 2 - dt 2 



with L^ds as the AdS characteristic length and the function f(z), 



(23) 



/W= (24) 

for which it holds that f(zo) = 0. In this coordinate system, the boundary is the z = 
slice and the geometry is capped off at z — z while keeping asymptotically AdS since 
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f(z) — > 1 as z — > 0. As it will be shown, different values of the parameters q\ and q 2 give 
rise to different IR behaviours for the holographic entanglement entropy and the two point 
functions. In this sense, the pure AdS space can be obtained by simply choosing qi = q 2 = 
in ([23]) while the simplest capped geometry, the non extremal soliton is defined by choosing 
qi = — g 2 = 1 in fflM . which results in f(z) = 1 — (j^ . 

In order to obtain the entanglement entropy of an static region A lying at the boundary 



of an asymptotically AdS spacetime, the Ryu and Takayanagi proposal [23| accomplishes 
the task throgh the computation of an static surface of minimal area in the bulk geometry. 
In AdS3, this means a static curve of minimal length 7, i.e a geodesic. The metric fT25|) on 
the fixed time slice is given by, 



ds 



T 2 

2 I _ n AdS 

SoiU=o - 



dz 



2 



+ f(z) dx 2 



(25) 



L/C*) 

Let us now consider the region A in the boundary given by an interval with length I 
along x direction. Choosing the origin in the center of this interval, the symmetry of the 
problem allows us to restrict to curves described by an even function z = z(x). Then, 
the length functional that we have to minimize to compute the holographic entanglement 
entropy reads, 



Length( 7 ) = 2 jf ^f- dx Jj^- + f(z) , (26) 

where z = dz/dx and z m satisfies that z(0) = z m while z(£/2) = e (see figure 7). Consider- 
ing the integrand of (1231) as a Lagrangian density C[z, z], one notices that it does not depend 
explicitly on x. Namely, the independence of C[z, z] on x leads to the conserved quantity 
H = p z z—C, where p z = dC/dz. In particular, one gets % = —(f(z)/z) (i 2 / / \z) + f '(z)) -1 / 2 . 
The constancy of Tl can be expressed by settling that z m = 1/H which allows to write, 

* = - ^ 2 , (27) 

as one notices that z < 0. Replacing (12"7|) in < HE\i , the geodesic length used in the formula 
for the computation of the holographic entanglement entropy (j20|) . acquires the form, 

f e z dz 

Length( 7 ) = -2L AdS / / ■ (28) 

J z m z y z m j yz) z 

In addition, the two point functions of scaling operators $ inserted at two spacelike 
separated points on the boundary of AdS are also considered. According to this correlator 
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, is obtained from the spacetime propagator <£jj olo between the corresponding points on 
the UV-cutoff boundary lying at z ~ e by simply writting, 

egW^eSo, (29) 

where T> is the distance between the inserted operators on the boundary and A is the 
scaling dimension of the operator $. The holographic propagator €jj olo , in the leading order 
of its semi-classical approximation, is given by a sum over geodesies connecting the 
points on the cut-off boundary, 



= e -mLength(Tu) ^ g -m Length^-p) j-gg^ 



Here, Length(7x>) is the length of the shortest geodesic connecting the two boundary 
points and m is the mass of the fields <fi dual to the scaling operatots $. In the AdS/CFT 
correspondence, the scaling dimensions of operators inserted at the boundary CFT are re- 
lated to the masses m of their dual fiels in the bulk geometry. For instance, when $ is an 

n 

scalar operator of large scaling dimension A, then it holds that A « mLAdS PI- 



A. Gapless limit: the pure AdS q\ = q 2 = 



The pure AdS space is obtained by choosing qi = q 2 = in equation ( 1241) . i.e, f(z) = 1. 
For convenience from here in advance, we simply write Lacls — L. Replacing f(z) = l into 
( |28l) and carrying out the integration one obtains, 



Length(7) = 2L 



log 



+ 7 



Z m ^ 2 



2L log 



2 



(31) 



As a result, the holographic entanglement entropy of a region A of length I reads (as one 
realizes that t = 2 z m , see figure 7 upper part left), 



S A 



Length(7) 



AG 



(3) 
N 



c _ r . 2z m c £ 

— 2L log = - log - 

dl e 3 e 



(32) 



Here, it has been used the celebrated relation between L/G$ and the central charge c 
of the boundary CFT given in [301 ] . 

3L 



2G 



(3) 
A' 



(33) 
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The equation (I3T[) also allows to compute the two point function of an scaling operator $ 
of large scaling dimension A inserted at two points separated by a distance T> = £. Making 
use of the equations ( )2"9"|) and (I3"0l) . this correlator amounts to, 

(jf =e -2A e -mLcngth( 7 ) = ( g 4) 

_ e -2A e -m2L\og£/e ^ 
~ e -2A e ~2A log£A = r 2A ? 

which is exactly the power-law decaying of the correlation functions expected for the 
ground state of a critical system which geometric dual is given by the pure AdS space. It 
is clear since 10] and the exposition above that, the behaviour of two point functions and 
entanglement entropy in pure AdS space are well represented in terms of an scale invariant 
MERA tensor network (see equations (fT5l) and (122|) in sections III Bl and II I II respectively) . 



B. Non extremal soliton q\ 



"<?2 



1 



The simplest AdS soliton is the non extremal soliton, obtained by taking qi = — q 2 
in equation (J2l|) . In this case, the function f(z) reads as, 

/(*) = \ (4 - z 2 ) ■ 

The computation of the geodesic length ([28]) for this geometry yields, 

a + V a 2 — z 2 



(35) 



Length(7) = 2L 



log- 



(36) 



where a 2 



In the limit z m <C zq, a ~ z m and ( 1361) reduces to ( 1311) . However, here 



we are interested in probing the entanglement entropy in the IR regions of the geometry 
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31 



32], i.e, when z m ^> z and a ~ Zq. In this limit, and requiring that the integral 



have a real value, the geodesic length amounts to, 



Length(7) = 2L 



log 



2L log 



2Zn 



(37) 



As a result, one finds a saturated length of 7 when probing the IR geometry of the 

translates into the well known saturation for the 



non extremal soliton, which through ([20 



entanglement entropy in gapped phases 23|, [33J , 



Sa = S ~ log 



2Zn 



(38) 
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In other words, in the non extremal soliton geometry, the entanglement entropy which 
corresponds to small I is similar to ( 152"]) and the geodesic has a turning point at z m . At 
variance, when the minimal curve probes the region z ~ z , its length saturates yielding 
the well known saturation of the entanglement entropy in gapped phases. Namely, when 
£ ^> zq, the minimal curve 7, after having probed all the length scales smaller than zq, 
finishes at two endpoints located at z = z (see figure 7, lower part, right) and the longest 
complete geodesic curve able to connect two boundary points is the 7 curve between two 
points separated at a distance V = 2 z with the turning point at z m = z Q . When V > 2 z , 
there are no geodesic curves directly connecting the boundary points, a signature of the 
vanishing of the two point functions for distances larger than 2 zq. 

The entangling nature of the IR fixed point exposed by the large I expansion of the 
minimal curve 7, can be interpreted in terms of the hybrid MPS-MERA tensor network 
states described in section III CI With this aim, a tensor network composed of a MERA 
curtain of tq = log £ layers, with £ = 2 zo is considered. On the top of the curtain, there 
is an MPS state \Q) with bond dimension A4 — 1 which corresponds to the trivial non 
entangled product state, 

|0> = |0> |0> -.- J0> . (39) 

Let us first to compute the entanglement entropy of a region A with i 3> £. In the MERA 
curtain, the causal cone of the region A is connected to the rest of the MERA sites through 
a number of traced out bonds that is proportional to log£ which is independent of £, i.e, 
Length(7 J 4) w log£. Using equation fT2~T]) one finds that the contribution to entanglement 
entropy given by the MERA curtain is, S^ ERA cx log£ ~ So that exactly matches the 
scaling of the entanglement entropy for the AdS non extremal soliton ( 13"5|) . On the other 
hand, the product state \Q) at the top of the network does not collect any contribution to 
the entanglement entropy of the region A since these contributions would be bounded by 
logA^ = logl = (see eq.()9])). As a result, the finishing of the geodesic curves 7 at two 
endpoints in z in the non extremal soliton, may be interpreted in terms of an MPS at the 
top of the MERA curtain which results to be the non entangled product state 

Next thing to consider is the correlation function of two scaling operators $(— E/2) and 
$(£/2) located at the bottom of the MERA curtain, £ e Tm = (tf |$(-£/2) $(^/2)|tf). After 
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To = log £ MERA contractions such as those described in section III Bl one obtains, 

(*|$(-^/2)$(^/2)|*)= (40) 
= (0| (Sm-£/2)) T0 (Smi/2)) T0 \Q) = 

= (^(nM-e/20H£/20\ty , 

where |\?) is the tensor decomposition of the system given by the hybrid MERA-MPS 
network. Noticing that A = — log//, and (jHJ), the correlator can be written as, 

(tf|$(-*/2)$(^/2)|tf}= (41) 
= T 2A (n\$(-£/2£)$(£/2Z)\n) = 

= r 2A £ c, = o , 

|A„|<1 

where in the last equality, it has been used that for the product state \ u = 0, V 
| A y | < 1. This implies that all the two point functions vanish for distances between the 
boundary points T> ^> £ = 2 zq which completely agrees with the geometric prediction using 
the AdS/CFT. Finally, it is worth to note here that the conformal behaviour obtained in 
section IIV Al is trivially recovered by the hybrid MERA-MPS tensor network when the size 
of region A is such that £ <C £. 



C. Extremal Soliton q\ > 1, qi > 1 



The extremal soliton geometry is defined by the class of harmonic functions 

f(z) = l-a(-)+/3 (- 



(42) 



where a = (91 + 92) and (3 = q\q<i- The geodesic length (|28|) in the AdS soliton geometry 
which corresponds to (J4"2]l reads as, 

dz 



Length(7) = — 2L z 



2L log 



Zm z^Zq - az z + gz 2 
2z — az + 2^/zq + £>,2 2 — ckZq, 



(43) 



with g = /3— (^j . In terms of the holographic prescription to compute the entanglement 
entropy (120]) . the hints of hidden topological orders encoded in the features of the long range 
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FIG. 7: Geodesies in AdS solitons. Upper half left: In Poincare coordinates, the geodesies con- 
necting two boundary points separated by distance t in the pure AdS3 spacetime are given by half 
circumferences with radius z m = £/2. Upper half right: An AdS soliton geometry is caracterized 
by a function f(z) that vanishes at a certain value zq of the radial coordinate z. This amounts 
to say that the holographic geometry given by equation (1231) ends at zq (discontinuous red line) 
so geodesies can no longer probe regions with z > zq. Lower half right: the geodesies in the non 
extremal soliton geometry end at two points lying at z = zq and Length^) ~ log2zo- Lower half 
left: the geodesies in the extremal soliton geometry, once ariving at z ~ zo, instead of finishing, 
bend-over the tip of the soliton geometry (red continuous line between points at z ~ zq) giving 
additional contributions to the geodesic length. 

entanglement, should be tracked down at the IR behaviour (the large i expansion) of the 
UV-independent piece of (T4"5I) . Hence, defining the variable K = z m /z and then expanding 
(f43|) for large fC one obtains, 




(44) 



2L log - + log 




where, 




(45) 
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£ = 2 zq and XJC < 1. It is importatant to note that A > for any selection holding 
simulaneously that q\ > 1, 92 > 1 and /(zo) = 0. 

In equation (jHJ) one finds two very different kind of terms: first there is an UV-dependent 
(e-dependent) logarithmic term that is responsible for the area law of the entanglement 
entropy when applying (1201) . In addition, there is also an IR term log(/C/ (1 — A/C)) which is 
quite a generic feature of holographic models with a gapped spectrum provided the growth 
of minimal surfaces at the IR region (see figure 7 lower part left). It is pertinent to note 
here that, if one reduces the size of the subsystem A (i.e, making I smaller), eventually, the 
curve 7 separates from the the tip of the soliton geometry at z = zq and the smaller size of 



A, now probes more directly, the UV region of the boundary field theory 32) . When this 
occurs, the IR term disappears. 

As in the case of the non extremal soliton, it is possible to interpret the entangling nature 
of the IR fixed point shown by the large K-t expansion of Length^), in the language of the 
hybrid MPS-MERA tensor network state described in section III CI Once again, a tensor 
network composed of a MERA curtain of r = log £ layers with £ = 2 Zq is considered. 
However, on the top of the MERA curtain, now lies an MPS state \i/jmps) with bond 
dimension M. > 1 that corresponds to a non trivial entangled state. 

A similar argument to the case of the non-extremal soliton for the entanglement entropy 
of a region A with £ ^> £, shows that the MERA curtain contribution to entanglement 
entropy is given by S A 4ERA oc log£ ~ Sq. Nevertheless, in this case, the state \iPmps) 
at the top of the network, provides an additional contribution S^ IPS to the entanglement 
entropy of A. This non zero contribution results bounded by \ogA4 (see eq.([9])) and the 
total entanglement entropy now amounts to, 

S A = SX IERA + S™ PS < log £ + 2 log M . (46) 

As the entanglement collection on the top layer now yiels a non- zero result, it is reasonable 
to interpret that the curve acting as the boundary of the causal cone of A, has an IR flat 
bend-over region near the MERA curtain top layer. This is also consistent with the growth 
of minimal surfaces at the capping region in holographic models with a gapped spectrum. 
It is thus precisely S^ PS , the term which may be identified with the IR contribution to the 
geodesic length in ( )44l) as far as it holds that, 

log ^<2 1ogM. (47) 

\-\K 
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The two point function of two scaling operators $(— 1/2) and $(£/2) located at the 
bottom of the tensor network, ^ns = (^|^ > ( — ^/2) $(£/2)|\I/) can also be derived from 
simple arguments. Once again, r = log£ MERA contractions are carried out in order to 
obtain that, 

(#|$(-£/2)$(£/2)|*) = (48) 

= ojmpsI (sm-e/2)r (smi/2)r \^ MPS ) = 
= ^) 2To (iJ MPS m-£/20H£m)\^MP S ) . 

Recalling that A = — log/x, and ([8]) one easily finds that, 

(*|$(-^/2)$(^/2)|*)= (49) 

= r 2A £ c ^ ( A ^)^ • 

|A„|<1 

Here it will be assumed that the sum in (j4"9~j) . may be approximated by the term with the 
highest \ u i.e, 

4ns = (*|*(-*/2)*(*/2)|¥) ^T 2A c 2 (A 2 )^ . (50) 

As in previous sections, it is interesting to compare this result with the two point function 
of two scaling operators of the boundary theory separated a distance £, when it is obtained 
through (jUJ) and the holographic prescription ( 1291 1301) . Proceeding in this way, one finds 
that, 

£W olog = e- 2A e- mLea z th ^ = (51) 

_ e -2A e -2A(logC/e+log/C/(l-A/C)) = 

2 A 



t- 2A 



1 - A/C 



By a direct comparison between ( 15 Up and (1511) . one realizes how similar are the holo- 
graphic structure of the two point functions offered from MERA-like tensor networks and 
the AdS / CFT. It is thus tempting to further identitfy some relevant parameters defining both 
approaches. Assuming that the order of the scaling dimension of operators <3? are 0(A) ~ 1, 
it is thus possible to identify the exponentially decaying terms and ((1 — A/C)//C) 2A 
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in ( 150]) and (l5~Tj) respectively, which yields, 



log- ^— — logA 2 = r?/C, (52) 

1 — A /C z 



with 



^ = -^!og^2- (53) 

As a result, equation flUj) may be expressed as, 

Length( 7 ) « 2L ^log | + 77 /C ^ . (54) 

The extensive term 77/C is consistent with the fact that, in general, any holographic 
geometry that is capped is a candidate for showing extensive behaviour provided the growth 
of the minimal surface 7 when it bends over the infrared region. The identification of 
the exponentially decaying terms posed in (1521) may be pushed even further in order to 
obtain a direct relationship between parameters of the tensor network and parameters of 
the holographic geometry. Indeed, considering the limit r))C <C 1 in (|52|) one finds that, 

V = ~\ lo S A 2 ~ ^p- + A , (55) 

and 

v\k^oo ~ Km !^p- + A = lim p + A = A . (56) 

As a result, one may conclude that the AdS/MERA duality provides a relationship be- 
tween the parameter A defining the extremal AdS soliton geometry and the eigenvalue A2 of 
the identity transfer matrix corresponding to the top MPS tensor. This connection can be 
written as, 

A~-ilogA 2 . (57) 



It is worthy to notice here that, all the identifications and conclusions raised from 
are only valid when A > 0. For the extremal soliton defined by qi — q% — 1, the parameter 
A is zero. Applying A = to ( 1541) one obtains that, 

Length( 7 ) « 2L log ~ . (58) 
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This result for the geodesic length corresponds to a a non-extremal soliton, which obvi- 
ously is a wrong conclusion because A = was applied where not allowed. Instead, when 
considering A = in equation (1H|) . one obtains, 

Length( 7 ) = 2L ^log | + log/C J = 2L log * , (59) 

which corresponds to the geodesic length in pure AdS3 or in AdS/MERA terms, to a 
pure scale invariant MERA tensor network. 



V. DISCUSSION AND CONCLUDING REMARKS 

By exploring some hints posed in jl^ . it has been explicitly shown that the entanglement 
entropy and the two point correlators in a type of hybrid tensor network state composed by 
a finite number of MERA layers and an MPS acting as a cap layer, resemble the behaviour 
of the holographic entanglement entropy and the holographic two point functions appearing 
in the AdS soliton geometries. Within the context of the AdS/CFT correspondence, the 
AdS soliton is a capped geometry representing a theory with a mass gap. It was thus 
expected that the entanglement at the IR point would be zero when considering models 
where no topological order is present. This is the case of the non extremal soliton, in which, 
by making use of the holographic formula for the entanglement entropy, no hints of long 
range entanglement were found. According to this, a tensor network representation of the 
non extremal soliton which is characterized by an MPS cap tensor that is a trivial product 
state has been proposed. At variance, for extremal soliton geometries, the holographic 
formula for the entanglement entropy leads to find nontrivial behaviours for the long range 
entanglement, which might be related to some underlying topological order In this case, 
the capping tensor in the dual tensor network results in a non trivially entangled matrix 
product state with finite, non zero bond dimension. 

Assuming the striking relationship between MERA and AdS/CFT, a connection between 
some parameters of the tensor network represenation and those parameters which define 
the holographic geometry of an AdS soliton is provided. Indeed, equation (|57|) relates the 
spectrum of the tranfer matrices of the cap MPS tensor with the parameters q\ and q-i 
defining the soliton in AdS. This kind of result may be useful in order to, via AdS/MERA, 
open the possibility to classify the gapped phases of ID-quantum many-body systems in 
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terms of their associated geometry. As posed in 34| . different topological phases of ID 
gapped systems with an MPS representation, can be characterized by the spectrum of the 
MPS tranfer matrices. 

The proposal in this paper could be extendable to higher dimensions. Namely, there are 
2D generalizations of MERA 35] and MPS (PEPS) [6] to reasonably construct 2D versions of 
the hybrid network proposed here. In 36J , authors presented the MERA tensor network for 
the Kitaev's toric code model [37], a topologically ordered phase in two spatial dimensions. 
The MERA tensor network for the toric code model is capped off by a top tensor which 
encodes the information about the topological order. Furthermore, and in a closer analogy to 
the tensor networks considered here, the toric code also possess a very simple representation 



in terms of a PEPS with a bond dimension M. = 2 [38[. Indeed, for any PEPS with a 
finite bond dimension, it is always simple to calculate the topological entanglement entropy 
defined in 



39]. 



Finally, in 40J], authors provide the gravity dual of a (2+1) dimensional field theory 
which in its IR limit flows to a pure Chern-Simons theory. Since the theory is gapped, the 
geometry is capped off at IR. To complete the model, D7-branes are located on the tip of 
the cigar-like soliton geometry, making the IR behavior different from that of the "naked" 
AdS soliton geometry. Concretly, from the D7-branes arises a Chern-Simons term, which 
gives rise to a topological entanglement entropy term through ffl9l) . It is thus tempting to 
apply the ideas presented here in the context of higher dimensional geometries i.e, finding 2D 
tensor network states with top tensors whose entanglement properties match those observed 
in the geometric duals of [40]. 
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